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Abstract 

A general expression for temperature-dependent magnetic susceptibility of 
quantum gases composed of particles possessing both charge and spin degrees 
of freedom has been obtained within the framework of the generalized random 
phase approximation. The conditions for the existence of dia-, para-, and ferro- 
magnetism have been analyzed in terms of a parameter involving single-particle 
charge and spin. The limit T — > retrieves the expressions for the Landau and 
the Pauli susceptibilities for an electron gas. It is found for a Bose gas that on 
decreasing the temperature, it passes either through a diamagnetic incomplete 
Meissner-effect regime or through a paramagnetic-ferromagnetic large magneti- 
zation fluctuation regime before going to the Meissner phase at T = Tbec- 
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1 Introduction 



We wish to draw attention to an interesting application of analytical expressions de- 
rived in Paper I [Tj for the particle-current and the transverse charge-current static 
susceptibilities of uniform gases of neutral quantum particles. We consider here M 
quantum particles, each carrying a charge (Ze) and a spin s, moving in a volume V 
in the presence of a neutralizing background of opposite charges and thus forming a 
quantum mono-plasma ("jellium"). The dimensionless temperature-dependent mag- 
netic susceptibility of such a fluid, which may be deduced as static limit from Ref.[2], 
Eq.(Al6)], is given by 



M-H uj^ 



H^o H q^o c q 
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where M, H, eo, c, Up, and xlil) denote, respectively, magnetization, magnetic field, 
permittivity of vacuum, speed of light, plasma frequency, and static transverse charge- 
current susceptibility. For the mono-plasma, one has ujp={Ze)'^n/{eom) with mass m 
and average number density n=N' /V of particles. 

Choosing a cartesian coordinate system {e^^, e^, e^} such that q=gej;, the static 
susceptibility of the transverse charge-current density is given by 

Xl{q) = J^^^,l clte-^'([jj,(t)t, JJ^(O)]) (2) 
^ ^0^' <! 1 + K^^:''Z'"' J n ) . (3) 
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where Jq ^ denotes the operator of transverse charge-current density. In the last line, 
Xx(q') has been approximated in terms of the transverse charge-current susceptibility 
X± (9) of ci fictitious similar system in which interactions between particles have 
been switched off 0. This generalized random phase approximation (GRPA) correctly 
accounts for the transverse electromagnetic shielding effects in systems of particles 
which carry a charge and/or a spin. The conventional RPA expression for the transverse 
'particle- current susceptibility of a system of neutral and spinless particles, Eq.(45) of 
Ref. pj, is easily recovered from Eq.([3]) (taking into account Eq.(jl]) below) for s=0 and 
Z ^ 0. 

There are three contributions to the operator of charge-current density [1] resulting 
in the transverse component 

Jj, = (Ze) Jq,. - ig7Sq,, - ^ E (e. • A^) 5iVq^k (4) 

with Jq denoting the number-current density, Sq= I]ko-o-' So-o-'ak^ak+qo-' the spin density, 
7=fl'/^B/^ the gyromagnetic ratio with g as the (^-factor and ii-Q=eh / {2me) the Bohr 
magneton, iVq the particle-number density and its fluctuation 5iVq=iVq — iVq ^, and 
Aq= / d^r e~"''"A(r) the vector potential associated with the moving charges. How- 
ever, in a uniform system only the first two terms on the right-hand side will have 
main contributions to the transverse charge-current susceptibility in Eq.([2]), because 
the contribution due to the vector potential Ak produced by the moving charges will be 
negligibly small in a non-relativistic gas. Moreover, cross correlations between Jq, ^ and 
Sq,j, will also vanish due to rotational symmetry. Therefore, we expect only two con- 
tributions to the static charge-current susceptibility: an 'orbital contribution' induced 
by the macroscopic current density due to moving charges and a 'spin contribution' in- 
duced by the magnetization- current density due to magnetic moments associated with 
particle's spin. 

Consequently, we express the transverse charge-current susceptibility of the nonin- 
teracting system, which is needed in Eq.Q, as 
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Here ^-iid Xyyl?) denote the ideal-gas susceptibilities of the particle-current den- 

sity and the spin density, respectively. The latter is easily shown to be proportional to 
the number-density susceptibility of an ideal quantum gas. The expressions for 
ideal-gas static susceptibilities x{.l) have been obtained in Paper I. For no- 

tational simplicity, superscripts on these quantities are being omitted now onwards. 



2 Magnetic Susceptibility 



The expression for the magnetic susceptibility within GRPA is finally obtained - using 
Eqs. dlD, (I3l) and (I5D - in the form 



where 



s(r) 



lim ^ 
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It is to be noted that S(T) may be interpreted as a static susceptibility, too, which 
describes the linear response of the magnetization M to the magnetic induction field 
B. From Eqs.([T]) and in conjunction with the general relation B=/io(l + Xm)H, 
where /io=l/(eoc^) is the free-space permeability, one finds 



, MB 2 



(8) 



In order to evaluate S(T), we read for the small-g region from Eqs. (34) and (36) of 
Paper I: 



q'^xiQ) 



and 



X±{(1) — 1 
n 
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Here P=l/{kBT) and A= exp[/3/i^(n, T)] denotes the fugacity with chemical potential 
firf{n,T), where 77= — 1, 0, +1 for the gas obeying FD, MB, BE statistics, respectively. 
Co= limAT^oo ^0/^) the fraction of particles occupying the zero-momentum single- 
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particle state, is given by Co=5r), 1 © (^bec ~ T) 

the unit-step function and Tbec the BEC critical temperature. Also, Ci,{x 
denotes the polylogarithm. 



with B(x) denoting 
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Inserting Eqs.([9]) and (fTOj) into Eq.Q, one gets 



sign(^ - 2^2) X cx) if Co > 



E{T) = \ "^'^''^^^^J^T '-o = u ^^^^ 



with the "Curie constant" 

To=^ L^ , 12 
12 eoc"^ kb 

and the parameter 

^ = s{s + l)g^ {m/m,f - (13) 

determining whether the gas will show paramagnetic > 0) or diamagnetic < 0) 
behaviour. 

Thus for neutral particles of nonzero spin {Z=0, s > 1/2 =^ ^ > 0, paramagnetic) 
one finds for the "Curie constant" the well known expression 

^ _ /iog(g + l)(c//iB)^n 

which does not depend on particle's mass m. The Curie law is recovered for a Boltz- 
mann gas since the ratio Ci/2(^A)/C3/2('7A) becomes unity. On the other hand, for 
spinless charged particles (Z > 0, s=0 =^ ^ < 0, diamagnetic), one finds 

which varies as m"^. The general expressions derived by us here give the diamag- 
netic behaviour of spinless non-degenerate charged bosons consistent with the result 
obtained in Ref. [5]. Furthermore, we find that a paramagnetic gas will show a fer- 
romagnetic ordering transition (xm=oo) at a Curie temperature Tc determined by 
S(Tc)=l as per Eq.Q. On the other hand, a diamagnetic Bose gas will exhibit the 
Meissner-Ochsenfeld effect (Xm= ~ !)> also called simply the Meissner effect, in the 
Bose-condensed phase (0 < T < Tbec) since |S(T)|=oo from Eq. (fTT]) . 

In Fig^H the (dimensionless) reduced susceptibility 

_^,^, sign(06„ 1 Ci(^A*) 

with A* = e'^^*^"'"^*^/^' is plotted as a function of the reduced temperature T*=k-QT/e^^ 
where eu=h^k^/ (2m) and ku=2 {67T'^n/(2s + 1))^^'^. It can be seen that S* for a BE (an 
MB) gas diverges at T=Tbec (^=0), while S*(0)=6 for an FD gas. 

Typical values of Tq* are of order 10"^. Hence the magnetic susceptibility of a Fermi 
gas will be given, to a good approximation, by 



"*(T^*) 
"*^2~'*) 



Xm = ' ^ sign(0 T* E*{T*) , (17) 



4 



3*(r) 




0.0 0.2 0.4 0.6 0.8 



Figure 1: Temperature-dependent reduced susceptibility for BE (full 
line), FD (dash-dotted), and MB (dashed) gas. The vertical gridline 
indicates the location of Tbec of the ideal Bose gas. 



because 1/Tq* ^ H*(T*) at all temperatures in contrast with MD and BE gases for 
which S*(T*) can take very large values. In order to demonstrate such a behaviour 
graphically, we have plotted function of temperature for the extraordinarily 

large value Tq*=10~^ and ^ < (diamagnetic behaviour) in Figj2] for all the three 
gases. Figure [3] demonstrates the behaviour of a paramagnetic gas (^ > 0) for the 
extraordinarily large value Tq*=10~^. For this value of Tq, whereas the FD gas remains 
paramagnetic at all T, the BE gas and the MB gas both show a ferromagnetic phase 
transition at a Tc which is just above Tbec for bosons and far below for MB gases. 
Also, the ferromagnetic Bose gas shows large magnetization fluctuations before going 
to the Meissner phase. 



If we consider the limiting behaviour T ^ for a Fermi gas, we find from Eqs. (fT6|) 
and ([nD together with Eq. ([12]) that 

^Ferrm^j. - 0) = sign(0 ^ (18) 

from which we obtain the celebrated results, see for example [^J, for an ideal electron 
gas: 

Y =--Y = - "'^^ (19) 

It is pertinent to mention here that Schafroth [7] had established before the advent 
of the Bardeen-Cooper-Schrieffer [8] theory of superconductivity that an ideal charged 
Bose gas (CBG) below Tbec shows the fundamental equilibrium manifestation of a 
superconductor, i.e. the Meissner effect. This implies that if the system is cooled 
from the high temperature (normal) phase, both the Bose-condensed and the super- 
conducting phases should appear simultaneously. The concept was supported further 
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Figure 2: Temperature dependence of susceptibility for diamagnetic 
systems: BE (full line), FD (dash-dotted), and MB (dashed) gas for 



rQ*=10 ^. Vertical gridline as in Fig. [H 
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Figure 3: Temperature dependence of susceptibility for paramagnetic 
systems: BE (full line), FD (dash-dotted), and MB (dashed) gas for 
Tq*=10~^. Vertical gridline as in Fig. [H 
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in works [9l [10] on the basis of gauge invariance together with the existence of the 
off-diagonal long-range order in systems which, in general, could be interacting. How- 
ever, the study of Koh [11] asserts that a CBG with short-range repulsion would show 
incomplete Meissner effect prior to onset of the BEG when complete Meissner effect 
would appear. It can be observed that our results depicted in Fig. [2] support such a 
conclusion. The behaviour in Xm^T plane at T > Tbec seems typically akin to that for 
type II superconductors in Xm^H plane wherein only partial expulsion of the magnetic 
flux takes place (i.e. \xm\ < 1) between the critical fields Hc^ and Hq^- Also, behaviour 
depicted in Fig. [3] suggests possibility of the role of the charged bosons possessing spin 
degree of freedom in understanding spin-triplet superconducting transition observed in 
ferromagnetic compounds, see, for example fi2\ [T3]. 



3 Conclusion 

A general expression for temperature-dependent magnetic susceptibility within GRPA 
has been obtained for the BE, FD, and MB gases by relating it to the long- wavelength 
properties of and x±{q) fo^' noninteracting particles. The expression involves both 
charge and spin of a particle. It is shown that the Landau (diamagnetic) and the 
Pauli (paramagnetic) susceptibilities are special cases (T ^ 0) of our more general 
expression. Also, the T-dependence of diamagnetic Xm for an ideal electron gas is in 
conformity with the numerical results obtained in [Hj although the comparison has not 
been depicted here. In summary, our work contains a unified approach for the three 
gases and provides a combined treatment for the study of the magnetic susceptibility 
for a system in which constituent particles possess both charge and spin. 
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